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Abstract 

Recent molecular and karyologic studies have significantly modified delimitation of Lilium. However, despite the importance 
of pollen evolution in the genus comprehensive studies with electron microscopy and evaluation of pollen evolution are 
lacking. Therefore, we studied pollen morphology in a sample of 65 individuals from 37 taxa covering all the sections 
distributed in the world, using scanning electron microscopy. Our collection of 49 individuals from 21 taxa covering all five 
sections in China was also included in the database. We found pollen tetrads in L. bakerianum. Based on present and 
previous studies, our results suggest that pollen from L. formosanum should be classified as a new type, Formosanum. 
Combined with morphological and molecular evidence, pollen sculpture patterns appear to reflect phylogenetic 
relationships and are useful for species or subsection delimitation. Based on a comprehensive survey and correlation with 
potential functional implications, we propose the following hypothesis: evolution of an exine sculpture shows pollen type 
trends from Martagon — » Callose — » Concolor — > Formosanum. The evolutionary trend regarding pollen sculpture and size 
could be related to selective pressure to adapt to environmental conditions. Pollen size and shape showed a significantly 
positive correlation with annual precipitation, and smaller pollen grains appear to adapt better in habitats with extreme 
conditions. Evolution trends in exine sculpture do not appear to be definitively correlated with pollen size and shape. 
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Introduction 

Approximately 1 1 0 to 115 Lilium species are distributed in the 
cold and temperate regions of the Northern Hemisphere [1,2], 
particularly in East Asia, the Himalayas and Hengduan Moun- 
tains, North America and Europe. A total of 55 species occur in 
China [2] . De Jong [3] and Patterson and Givnish [4] consider 
southwest China and the Himalayas to be the center of origin of 
this genus. 

Classification of this genus has been historically complicated. 
Several classifications for Lilium have been proposed based on 
morphological characters. Detailed studies have been performed 
for East Asia, European and North American Lilium species [5-7] . 
Based on 1 3 morphological characters and two germination types, 
Comber divided this genus into the following seven sections: 
Martagon Rchb., Pseudolirium Endl. which is limited to North 
America, Liriotypus Asch. and Graeb. which is distributed across 
Europe and the Caucasus, Archelirion Baker, Sinomartagon Comber, 
Leucolirion Wilson, and Daurolirion Comber, representing the most 
widely accepted taxonomical divisions [8]. Wang and Tang 
recognized sect. Lophophorum (Bur. et Franch.) Wang et Tang out of 
sect. Sinomartagon Comber and included campaniform-flowered 
species [9]. Liang [10] and Haw [11] modified sect. Lophophorum to 
accommodate the Nomocharis-Yikc Lilium species in sect. Sinomartagon 



Comber. Chinese species were divided into five sections: Martagon, 
Archelirion, Sinomartagon, Leucolirion and Lophophorum. 

Recently, molecular phylogenetic analyses and chromosome 
techniques have improved the understanding of several groups 
within the genus and modified the phylogenetic position of 
Comber's classification, such as placement of sect. Daurolirion 
Comber in sect. Sinomartagon, L. henryi in subsect. Leucolirion 6b, 
modification of sect. Lophophorum and relationship confirmation in 
sect. Liriotypus [12-13,16-21]. Preliminary research found that 
sect. Sinomartagon, which mainly occurred in China, was compli- 
cated and polyphyletic. As indicated by Patterson and Givnish [4], 
intercontinental dispersal details of the genus Lilium are not yet 
clear. The division of subsect. Sinomartagon 5c and sect. Lophophorum 
are controversial and will require further research. Nishikawa 
et al. [18,19] suggested that L. henryi be classified into subsect. 
Leucolirion 6a and that it showed similar morphological features 
with L. rosthornii, thereby demonstrating that the phylogenetic 
position of L. rosthornii needs further study. 

There are few relevant studies regarding pollen morphology 
which defines the taxonomic and reflect the evolution of the genus. 
According to the description within Lilium by Baranova [22] based 
on the number, shape and arrangement of columellae that form 
the muri, there are three morphological types of pollen: (1) 
Martagon (muri formed by rectangular columellae); (2) Callose 
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(muri formed by rounded columellae); and (3) Concolor (muri 
formed by separated rounded and polygonal columellae). Previous 
studies found that most Lilium species have single pollen grains. 
However, pollen tetrads were found in L. sempervivouleum H. Lev. 
and L. amoenum E. H. Wilson ex Sealy, and the size and sculptural 
elements confirmed the taxa as two subspecies in L. sempervivoideum 
[23]. The pollen morphology of L. lophophorum (Bur. et Franch.) 
Franch., L. henrici Franch., L. souliei (Franch.) Sealy and L. nanum 
Klotz. et Garcke supported placement in Lilium, which differs from 
Nomocharis in apture and sculptural elements, and showed an 
evolutionary aperture trend from monocolpate to porate [24]. 
Muratovic et al. [25] showed that two related European species, L. 
bosniacum and L. carniolicum, share similar pollen morphology. In 
addition, pollen morphology of some Chinese species under 
scanning electron microscopy (SEM), including 9 species described 
by Li and Qin [26], 10 species and 3 cultivars by Zhang et al. [27], 
and 12 species and 6 cultivars by Wu et al. [28], could provide 
taxonomic implications within Lilium:pollen has not only the 
commonness of genus, but also the specificity of single species. 
Interspecific pollen size and morphological characteristics has 
some difference, which has a certain reference value for the 
classification of Lilium. For example, there are differences in pollen 
size among L. cernnum, L. lancijlium and L. pumilum. Also there are 
significant differences in pollen ornamentation and morphology 
between L. leucanthum and its variety L. leucanthum var. centifolium 
from Qinling Mountains. Results obtained by Wang et al. [29] 
and Liu [30] indicate that pollen size parameters of L. pumilum and 
L. concolor from different provenances showed different degrees of 
variation. Determining whether this is a universal phenomenon in 
other species or if it is an important characteristic of species from 
different provenances requires a large amount of data. Therefore, 
a comprehensive study of Lilium pollen morphology in China, 
particularly a comparative study between populations or prove- 
nances using SEM, was needed. Environmental constraints act on 
patterns of differentiation, as all character states are the 
consequence of interactions between phylogenetic and environ- 
mental constraints [31]. Pollen could be considered a functional 
unit, with the exine structure as a compromise between four main 
functions [32]: protective, harmomegathic (the ability to absorb 
bending stresses, as occur during desiccation), reservoir (the role in 
producing an adhesive surface or as recognition substances) and 
clustering. Based on a comparative survey in the Rhapontkum group 
(Asteraceae, Centaureinae), Hidalgo et al. [33] hypothesized that 
occurrence in habitats with extreme conditions appears to be the 
main factor in the evolution from one form of pollen to another. 

By combining previous and present results, our aim is to provide 
palynology information regarding the extent of pollen morpho- 
logical diversity within and between species/ sections, link this 
information with current phylogenetic data, and contribute to a 
comprehensive understanding of the systematics and evolution of 
Lilium by: (1) verifying pollen types in Chinese species; (2) 
discussing these findings within a phylogenetic framework; and 
(3) correlating pollen morphology and size with habitat. 

Materials and Methods 

Ethics Statement 

For conservation reasons, we collected only small amounts of 
Lilium pollen. Our activities did not have any adverse effects of 
Lilium populations, and permission for collecting was obtained 
from National Engineering Research Center for Flowers of China 
(No. 35 Tsinghua East Road Haidian District, Beijing, P.R. 
China). 



Table 1 shows the origin of all species studied which were 
collected fresh. Our observations for 21 taxa represented by 49 
samples were also included in the database. Division into the 
sections was based on studies by Comber [8], Gao et al. [13], 
Nishikawa et al. [18,19] and Du et al. [16]) as follows: two species 
and one variety of sect. Martagon, nine species and three varieties of 
sect. Sinomartagon including L. dauricum Ker-Gawl (sect. Daurolirion 
Comber), eight species and one variety of sect. Leucolirion including 
L. brownii F. E. Brown ex Miellez (sect. Archelirion Comber), L. henryi 
Baker and L. rosthornii Diels (sect. Sinomartagon Comber), one species 
of sect. Archelirion, two species of sect. Liriotypus, three species of 
sect. Pseudolirium, and six species and one variety of sect. 
Lophophorum (Table 1). 

To compare pollen features at the genus level, a database with 
results of the present study and those of Halbritter [34], Kosenko 
[35], Liang et al. [23,24], Muratovic et al. [25], Wu et al. [28], 
Zhang et al. [27] and Zhang et al. [36] was established (Table 1). 
Because the results do not permit a calculation of all the 
parameters cited above, our observations from 49 samples 
representing 2 1 taxa was used only for statistical analyses. 

Morphological observations were conducted using SEM. Pollen 
was directly attached to double-sided adhesive tape and examined 
under the microscope to locate the pollen. The sample was then 
taped to the object stage. Following gold spray coating, 
observation and image acquisition were conducted using a Hitachi 
S-3400 scanning electron microscope following the method 
described by Avetissian [37]. All microscopy procedures were 
performed at the Biotechnology Centre, Beijing Forestry Univer- 
sity, China. 

Biometric measurements were made using Image-Pro Plus 6.0 
(Media Cybernetics, USA). For each sample, 10 (5 for L. brownii 
from Shennongjia, Hubei Province) fully developed pollen grains 
were measured. Parameters considered were the polar axis (Pj, 
equatorial diameter (E), PIE ratio, lumina diameter and muri 
width (Table 1). The description of pollen morphology was based 
on the shape and sculpturing classifications of Baranova [22]. 
Arithmetic mean and standard deviations are shown for each 
parameter. 

To investigate pollen evolution, pollen features from the 
database were superimposed onto a phylogenetic framework 
constructed by internal transcribed spacer (ITS) sequences. We 
selected 23 taxa, including Cardiocrinum gigantum as an outgroup 
based on preliminary study [16]. The other 16 taxa, including an 
outgroup of Northolirion bulbiliferum, were cited [13,18-20]. All of 
the sequence data from the ITS regions were double-checked 
visually, edited in BioEdit 5.0.6 [38], and aligned in Clustal X 1.83 
[39] using default settings. Phylogenetic analyses were performed 
using maximum parsimony (MP) and maximum likelihood (ML). 
The MP tree was constructed with PAUP* 4.0b 10 [40]. ML 
phylogenetic analysis was performed using RAxML v. 7.0.4 with 
unique model parameters for ITS sequences [41]. A general time 
reversible model (GTR) was applied with discrete gamma 
distribution. Bootstrap pseudo replicates were performed 1000 
times using the fast bootstrapping option and the best scoring ML 
tree. Phylogenetic trees were visualized using Treeview [42] . The 
best scoring tree was visualized with FigTree vl.3.1 (http://tree. 
bio.ed.ac.uk/). 

SPSS 18.0 was used for statistical analyses. One-way analysis of 
variance (ANOVA) was used to evaluate whether or not 
differences in pollen characteristics from different provenances 
were significant. In those cases in which ANOVA revealed 
significant differences, a least significant difference (LSD) test was 
performed (Table 2). Correlation analysis was used to evaluate the 
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relationship between pollen parameters and habitat (Our collec- 
tion of 21 taxa, including 49 samples) (Table 3 and 4). 

Results 

Numerical results of the palynological study are summarized in 
Table 1. According to the Lilium classification by Comber [8], 
Nishikawa et al. [18,19] and Du et al. (2013) [16], palynological 
results for each section are as follows: 

Martagon 

Representatives of Martagon (Figure 1A— F) have pollen size 
ranging from L. distichum Nakai et Kamibayashi (91.56x34.57 u,m, 
P/E of 2.66) to L. tsingtauense Gilg (104.70x44.52 urn, P/E of 
2.36). The shape can be long-ellipsoidal. L. tsingtauense and L. 
distichum pollen show similar exine sculptures with a slight 
difference in that the intersection of the muri presents a 
protuberant and enlarged columellae in L. tsingtauense pollen 
(Figure 1E-F). L. tsingtauense and L. distichum correspond to 
Martagon type pollen, which is in agreement with Baranova [22]. 

Sinomartagon 

Triploid L. lancifolium pollen is relatively larger than the other 
diploid species in section Sinomartagon. The shape can be ellipsoid 
to long-ellipsoidal. Species belonging to Sinomartagon (Figure 1M-P; 
2A-P; 3A-P; 4A-P) show a decreasing tendency in size with 
regard to pollen features from subsect. 5a to subsect. 5c as follows: 
L. lancifolium Thunb (111.20x41.16 |jm, P/E of 2.72); L. davidii 
Duchartre ex Elwes (93.09x33.55 |xm, P/E of 2.79); L. davidii var. 
willmottiae (E. H. Wilson) Raffill (91.89x33.59, P/E of 2.75); L. 
duchartrei Franch (85.07x33.31, P/E of 2.56); L. leichtlinii Hook. f. 
var. maximowiczii (Regel) Baker (82.81x33.04, P/E of 2.52); L. 
pumilum DC (87.39x32.66, P/E of 2.67); L. concolor Salisb. var. 
pulchellum (Fisch.) Regel (70.77 x30.58, P/E of 2.33); L. wardii Stapf 
ex Steam (87.77x35.26, P/E of 2.49); and L. taliense Franch 
(74.70x29.70, P/E of 2.50) (Table 1). 

SEM of pollen from this section showed differences in exine 
sculpture. L. concolor var. pulchellum clearly shows Concolor type 
pollen (Figure 4E-H), whereas pollen from other Sinomartagon 
species show rectangular columellae that correspond to Martagon 
type. 

Leucolirion 

Representatives of Leucolirion (Figure 1G-L; 5C-P; 6A-P; 7A-B) 
have pollen sizes ranging from 87.56x38.69, P/E of 2.28 for L. 
sulphureum Baker apud Hook. f. to 120.13x53.16, P/E of 2.29 for 
L. brownii from Shennongjia in the Hubei Province. The exine is 
ornamented with reticulation, and the shape is ellipsoidal to long- 
ellipsoidal. Pollen from L. brownii, L. henryi and L. rosthornii 
correspond to the Martagon type (Figure 1G-L; 5C^J). The 
pollen exine of L. regale Wilson, L. leucanthum (Baker) Baker, L. 
sargentiae Wilson and L. sulphureum show rounded columellae, 
corresponding to the Callose type (Figure 5M-P; 6A-P). Pollen of 
L. formosanum Wallace does not correspond to any of the pollen 
types proposed by Baranova (1985), and the exine surface is 
reticulate with solid muri or irregular-rugulate (Figure 7A-B). We 
suggest creating a new category of Formosanum type pollen. 

Lophophorum 

The division of Lophophorum remains controversial. Due to 
resource constraints, we selected L. bakerianum as the only 
representative. In a preliminary study, L. bakerianum was placed 
in the subsect. Lophophorum III [16]. Pollen tetrads were found in L. 
bakerianum (Figure 5A-B), observed as blunt quadrangles in a polar 



view with a pollen size of 78.40x79.71 um, P/E of 1.03. The 
exine ornamentation was reticulate with muri formed by separated 
rounded columellae, corresponding to the Callose type. This is the 
first study to report that pollen grains of L. bakerianum axe tetrads. 

Discussion 

Pollen Morphology and Phylogenetic Implications Based 
on ITS Sequences 

Within the genus Lilium, several diversified pollen exine 
sculptures were observed. Pollen sculpture patterns appear to 
reflect phylogenetic relationships and are useful for species or 
subsection delimitations within a section. 

The phylogenetic tree was resolved into six groups based on ITS 
sequences (Figure 8). Group I is comprised of three clades. Present 
and previous molecular analyses [18,19] and cytology studies 
[21,43-44] indicated that sect. Sinomartagon is complicated and 
polyphyletic, as supported by the pollen data (Figure 8). Species of 
sect. Sinomartagon (5a&5b) and L. formosanum and L. brownii (sect. 
Leucolirion 6b) formed a clade with strong support ([ML] BS = 81, 
[MP] BS = 77). Sect. Sinomartagon species show two pollen types, 
Martagon and Concolor. Differences in the shape and arrange- 
ment of the types of columellae could have implications on species 
delimitation. L. concolor var. pulchellum clearly shows Concolor type 
pollen (Figure 4G-J), and the remainder show Martagon type 
pollen. Meanwhile, three subtypes were showed within Martagon 
type pollen (Figure 8A, B, C). Exine sculptures of L. davidii and L. 
davidii var. willmottiae showed spherical but dorsally and ventrally 
compressed, and loosely arranged columellae with protuberant 
and enlarged columellae at the intersection of the muri (Figure 3B- 
N; 99B). Muri of L. leichtlinii var. maximowiczii have regular long- 
rectangular columellae (Figure 30-P; 8C). 

In addition, L. formosanum and L. brownii [Leucolirion 6b) formed a 
strongly supported clade ([ML] BS = 87, [MP] BS = 86) (Figure 8). 
In the cytological study, both L. formosanum and L. brownii lacked 
intercalary satellites on the first two pairs of chromosomes [43,46] 
However, the pollen of L. formosanum corresponds to the new 
Formosanum type, while L. brownii clearly shows Martagon type 
pollen. 

In our analysis, L. taliense and L. wardii ([ML] BS = 80, [MP] 
BS = 85) formed a clade with L. amoenum, L. sempervivoideum and L. 
bakerianum ([ML] BS = 79, [MP] BS = 84) that was moderately 
supported ([ML] BS = 56, [MP] BS<50) (Figure 8). In cytological 
studies, the L. taliense karyotype resembles that of L. wardii [21,43]. 
Morphologically, L. taliense and L. wardii showed a close 
relationship with each other, with a black line in the central 
groove of the tepals [2]. Pollen evidence from our study supports 
this result, as the pollen morphology of the two species is generally 
similar Martagon type (Figure 4K— P). Morphologically, L. 
amoenum, L. sempervivoideum and L. bakerianum show nodding or 
horizontal campanulate flowers, nectaries of inner tepals neither 
papillose nor with fimbriate projections, tepals without blotch at 
the base adaxially, and a papillose stem [2] . The L. bakerianum var. 
delavayi karyotype resembles that of L. sempervivoideum [43] . Pollen 
tetrads are present in L. amoenum, L. sempervivoideum [23] and L. 
bakerianum (Figure 5A-B). The similarity in pollen morphology also 
suggests a close relationship with this clade. 

The Martagon clade was monophyletic with strong support 
([ML] BS = 98,[MP] BS = 96), which is consistent with previous 
studies by Gao et al. [13] and Nishikawa et al. [18,19]. 
Morphologically, they all show whorled leaves. Although the 
Martagon clade shows Martagon type pollen, there are minor 
differences among species in the Martagon clade in the protuberant 
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Table 4. Correlation coefficient between the main characteristic of pollen morphology and the geographical and climate factors 
of provenances. 



Traits 


Latitude (°) 


Longitude (°) Altitude (m 


Annual average temperature (°C)Yearly precipitation (mm] 


Annual average sunshine (h) 


P 


0.046 


0.082 


-0.013 


0.172 


0.413* 


-0.013 


E 


-0.105 


-0.103 


0.041 


0.014 


-0.129 


0.249 


P/E 


0.133 


0.178 


-0.137 


0.160 


0.442** 


-0.251 


Lumina 


0.010 


-0.073 


0.039 


0.073 


0.024 


-0.091 


Muri 


0.126 


0.048 


-0.035 


-0.199 


-0.189 


0.134 



Note: * represents the significant difference (p<0.05), and ** represents the extremely significant difference (p<0.01}. 
doi:1 0.1 371 /journal.pone.0087841 .t004 



and enlarged columellae at the intersection of the muri in L. 
tsingtauense (Figure 1E-F, 8K). 

Group II was weakly supported ([ML] BS<50, [MP] BS<50) 
and consisted of two Lophophorum clades. Lophophorum clade I 
consisted of L. lophophorum and L. nanum ([ML] BS — 77 and [MP] 
BS — 7 1). In palynological research, L. lophophorum and L. nanum are 
monocolpate or 2-3 porate and show similar pollen morphologies 
(Figure 8L) [24]. Representatives of L. henrkii in subsect. 
Lophophorum II show Concolor type pollen with much smaller 
rectangular columellae between two separate irregular round 
columellae compared with L. concolor var. pulchellum (Figure 8M) 
[24]. Morphologically, representatives of the Lophophorum clade I 
show nectaries on the inner tepals with fimbriate/cristate 
projections on both surfaces, while representatives of subsect. 
Lophophorum II do not have this feature. The distinct pollen 
morphology of the species suggests a distinct relationship to 
support the molecular data [2]. 

In Group III, the Pseudolirium clade was resolved with strong 
support ([ML] BS = 77, [MP] BS = 86). PoUen morphologies from 
representatives of sect. Pseudolirium were similar, with loosely 
arranged columellae and protuberant and enlarged columellae at 
the intersection of the muri (Figure 8N) [35]. 

Group IV, consisting of L. kesselringianum, L. bosniacum and L. 
camiolicum, was robustly supported ([ML] BS = 97, [MP] BS = 97). 
The three representatives of European species (sect. Liriotypus) 
showed Martagon type pollen with densely arranged columellae 
forming the muri (Figure 80) [25]. 

Group V was only comprised of L. duchartrei [Sinomartagon clade I) 
(Figure 8). Based on morphology, Comber classified it in 
subsection Sinomartagon 5a. However, in our ITS analyses, L. 
duchartrei showed a distant relationship with subsect. 5a species as 
well as other sect. Sinomartagon species (Figure 8), which is in 
accordance with previous studies [13,18-19]. L. duchartrei {Sino- 
martagon clade I) karyotypes have been shown to be dissimilar to 
subsect. 5a Sinomartagon clade I species [21,43]. L. duchartrei pollen 
resembles that of L. taliense and L. wardii in exine sculpture 
(Figure 4C-D; K-P). Given that the phylogenetic position of L. 
duchartrei is uncertain, classification as a subsection of sect. 
Sinomartagon may be reasonable. 

In Group VI, L. speciosum var. gloriosoides in sect. Archelirion 
formed a weakly supported group consisting of representatives of 
Leucolirion 6a ([ML] BS<50 and [MP] BS<50). Leucolirion 6a 
comprised two subclades ([ML] BS = 74 and [MP] BS = 75). L. 
speciosum var. gloriosoides showed a distinct relationship with L. henryi 
and L. rosthornii in the ITS tree, although they all have Martagon- 
type pollen (Figure 8). In addition, L. henryi and L. rosthornii, which 
have orange reflexed flowers with prominent papillae and 
pubescent nectaries, were resolved with strong support ([ML] 
BS = 98, [MP] BS = 100). Species of L. leucanthum, L. sargentiae, L. 



sulphureum and L. regale formed a clade with strong support ([ML] 
BS = 95, [MP] BS = 95). Pollen morphology from our study 
supports this result, as the pollen morphology of L. henryi resembles 
L. rosthornii as Martagon type (Figure 5C-H), while L. leucanthum, L. 
sargentiae, L. sulphureum and L. regale have all generally similar 
Callose-type pollen (Figure 5K-L; 6A-P). 

Within Lilium, we noted that one section may show two or more 
different pollen types (Figure 8). For instance, sect. Sinomartagon 
shows two exine sculpture pollen types, Martagon and Concolor, 
and the former contains three different subtypes. No unique and 
uniform characters can be used to clearly distinguish sections or 
some subsections due to shared or overlapping characters. 

Trends in Exine Sculpture Evolution and Factors 
Implicated in Pollen Evolutionary Trends 

Changes in pollen morphology occur in response to selective 
pressure. Pollen grains should be considered a functional unit, with 
the exine structure as a compromise between the following four 
main functions [32]: protective, harmomegathic, reservoir and 
clustering. These functions are in response to physical components 
of the environment, such as water, nutrients, temperatures and 
growing season. 

The genus Lilium is distributed in East Asia, the Himalayas and 
Hengduan Mountains, Europe and North America [1,2]. 
Patterson and Givnish [4] concluded that Lilium evolved in the 
Himalayas and then dispersed into Eurasia and North America, 
although intercontinental dispersal details are not clear [4] . In the 
karyotype study, the relatively primitive group (sect. Lophophorum 
(Bur. et Franch.) Wang et Tang) have no intercalary satellites on 
the first two pairs of chromosomes. As expected, the genus 
JVotholirion which is commonly considered to be closely related to 
and more ancient than Lilium have also been found to lack 
intercalary satellites [4,21,43,45,46], these species are distributed 
mainly in the Himalayas, which makes this region more likely as 
the origin of the genus Lilium. Besides, Stewart [45] pointed out 
that secondary constrictions correlated with chromatin distribution 
as well as with geographic distribution. In his study, most East 
Asian and Himalayan species were found to have intercalary 
satellites on the first two pairs of submedian centromeres, while all 
North American species and European species lack such feature 
[21,25,43,45,]. Hence, the situation among the Himalayan and 
Henduan Mountains species appears more complicated since they 
seem to have both types, even within the same section indicating 
that the Hengduan Mountains are more likely as the differenti- 
ation centre of the genus [43,45,47]. Pollen evidence from our 
study supports this hypothesis, as the pollen morphology of the 
species from the Himalayas and Hengduan Mountains contain all 
of the pollen types in the genus Lilium. It should be emphasized 
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Figure 1. SEM photographs of pollen grains of Lilium. Figs A, B, L. distichum (Jia 07001). Figs C, D, L distichum (Du 11005). Figs E, F, L. 
tsingtauense. Figs G, H, L. brownii (Wang 1 1 006). Figs I, J, L. brownii (Du 09004). Figs K, L, L. brownii (Du 11014). Figs M, N, L lancifolium (Du 08001 ). Figs 
0, P, L. lancifolium (Wang 11001). Scale bars: 10 urn. 
doi:1 0.1 371 /journal.pone.0087841 .g001 



that the Formosanum type of pollen can be found in the genus 
JVomocharis, which is restricted to the Himalayas and Hengduan 
Mountains and it has been suggested that it be included in Lilium 
[13—16]. In contrast, the other regions possess a relatively singular 



pollen type, which led us to investigate Lilium pollen-type 
evolution. 

Based on a pollen study of 30 genera in Liliaceae, including 69 
species, using light microscopy (LM), Nair and Sharma [46] 
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Figure 2. SEM photographs of pollen grains of Lilium. Figs A, B, L lancifolium (Du 09001). Figs C, D, L. lancifolium(Du 11001). Figs E, F, L 
lancifolium (Du 09002). Figs G, H, L lancifolium (Du 11002). Figs I, J, L, lancifolium (Du 11003). Figs K, L, L. lancifolium (Du 09003). Figs M, N, L. 
lancifolium (Wang 11002). Figs 0, P, L. lancifolium (Du 11004). Scale bars: 10 Jim. 
doi:1 0.1 371 /journal. pone.0087841 .g002 
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Figure 3. SEM photographs of pollen grains of Lilium. Figs A, B, L. dauricum. Figs C, D, L davidii (Wang 1 1002). Figs E, F. L davidii (Du 1 1008). 
Figs G, H, L. davidii (Wang 1 1003). Figs I, J, L. davidii (Du 1 1009). Figs K, L, L. davidii var. willmottiae (Jia 07002). Figs M, N, L. davidii var. willmottiae 
(Wang 1 1004). Figs O, P, L leichtlinii var. maximowiczii (Du 12001). Scale bars: 10 |im. 
doi:1 0.1 371 /journal.pone.0087841 .g003 
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Figure 4. SEM photographs of pollen grains of Lilium. Figs A, B, L leichtlinii var. maximowiczii (Du 11010). Figs C, D, Lduchartrei. Figs E, F, 
Lpumilum. Figs G, H, Lconcolor var. pulchellum (Du 1 101 1). Figs I, J, L.concolor var. pulchellum (Du 1 1012). Figs K, L, L taliense (Jia 09002). Figs M, N, 
L.taliense (Du 12002). Figs 0, P, Lwardii. Scale bars: 10 \im. 
doi:1 0.1 371 /journal. pone.0087841 .g004 
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Figure 5. SEM photographs of pollen grains of Lilium. Figs A, B, Lbakerianum. Figs C, D, Lhenryi (Du 1 1015). Figs E, F, L.henryi (Du 12004). Figs 
G, H, Lrosthornii (Jia 1 0001 ). Figs I, J, Lrosthornii (Du 1 2005). Figs K, L, L. regale (Du 1 1 01 5). Figs M, N, L. regale (Wang 1 1 007). Figs 0, P, L. regale (Du 
11016). Scale bars: 10 \im. 
doi:10.1371/journal.pone.0087841.g005 
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Figure 6. SEM photographs of pollen grains of Lilium. Figs A, B, L regale (Du 12006). Figs C, D, Lsargentiae (Jia 1 1001). Figs E, F, Lsargentiae (Du 
1 1017). Figs G, H, Lsargentiae (Du 1 1018). Figs I, J, Lsargentiae (Wang 1 1008). Figs K, L, L. sargentiae (Wang 1 1009). Figs M, N, L, leucanthum. Figs 0, P. 
L. sulphureum. Scale bars: 10 |im. 
doi:1 0.1 371 /journal.pone.0087841 .g006 

proposed the scheme of pollen exine sculpture evolution (Liliaceae) Regarding pollen morphology of the genus Momochmis and Lilium 
as follows: ornate — * retipilate — > reticulate — » ring-shaped under LM and SEM, Liang and Zhang [24] suggested the 
reticulate — > scattered reticulate — » exineless ornamentation. following evolution trend of pollen exine sculpture in the genus 
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Figure 7. SEM photographs of pollen grains of Lilium. Figs A, B, L. formosanum. Scale bars: 10 |j.m. 
doi:1 0.1 371 /journal.pone.0087841 .g007 



Lilium: ornate — > retipilate — » reticulate, which is consistent with 
the hypothesis of Nair and Sharma [46] . In addition, retipilate was 
divided into three exine sculptures to form three Lilium pollen 
types, including Martagon, Callose and Concolor. These results 
led us to propose the following hypothesis of pollen evolution in 
exine sculptures within Lilium: Martagon type — » Callose type —> 
Concolor type — » Formosanum type. This also shows a reduction 
and simplification tendency. 

In monocots, reduction and simplification of the exine structure 
is even more extreme than in dicots and culminate in certain 
evolutionary orders [48]. Zavada proposed major evolutionary 
trends of wall structure types in monocots as follows: primitive 
tectate-columellate (perforate or imperforate) wall structures give 
rise to monocotyledonous atectate or granular walls and eventual 
extreme reduction of the exine, which may be completely absent 
[48]. All morphological characteristics are the consequence of 
interactions between phylogenetic and environmental constraints 
[32]. There is no doubt that environmental constraints act on 
patterns of differentiation in Lilium [49] . 

In addition to environmental factors, are there other factors 
such as pollinator affecting pollen differentiation? Numerous 
recent findings on the intricate flower pollinator networks suggest 
that differentiation of flora structures and functions in plants are 
tightly connected with those of pollinators [50-53]. Flowers of 
Lilium are insect-pollinated. Brantjes and Bos's [54] report that 
both diurnal and nocturnal hawkmoths visit L. martogon. Skinner 
[55] observed pollinators of 1 3 American species of Lilium and 
recognized five types of pollination: butterfly pollination (L. 
humboldtii, L. kelleyanum and L. kelloggii), butterfly and hummingbird 
pollination (L. wigginsii, L. pardalinum and L. pardalinum ssp. vollmeri), 
hummingbird pollination (L. columbianum, L. occidentale and L. 
bolanderi), hummingbird and bumblebee pollination (L. parvum and 
L. maritimum) and hawkmoth pollination (L. washingtonianum and L. 
parryi). Additional studies reported various flower visitors, includ- 
ing fritillaries to L. concolor var. pulchellum [56], swallowtails to L. 
dauricum [56], nocturnal hawkmoths to L. formosanum [57] and L. 
japonicum [58], and hawkmoths and swallowtails to L. auratum 
[56,59]. Nevertheless, wild Lilium species are self-compatible. 
Autogamy might have ensured reproductive success for the species 



in the environments where pollinators are comparatively rare [60] 
and the growing season short as well as where other conditions are 
selective. Pollen of L. martagon and L. dauricum with different 
pollinators are Martagon type (Table 1) [54,56]. Thus, for Lilium, 
the pollinators have little effect on changes of pollen morphology 
and size. Moreover, in Sonneratiaceae, the genus Duabanga has 
much smaller pollen grains than Sonneratia while both genera are 
bat-pollinated, hence absence of correlation with pollinator [32]. 
However, how floral traits as well as pollen grains of lilies are 
adapted to various pollinators remains to be studied in more detail. 

As mentioned, environmental factors constitute a selective 
pressure that may produce changes in the pollen structure. The 
number and arrangement of the types of columellae could 
influence exine flexibility. Exines with loose arrangement of round 
columellae (Callose, Concolor) are relatively more flexible 
compared to exines with tight, rectangular columellae (Martagon). 
This exine structure corresponds to harmomegathic function and 
the ability to absorb bending stresses that may occur during 
desiccation [32]. 

Relationship between Pollen Size and Conditions in the 
Genus Lilium 

Species at high elevations with relatively extreme conditions in 
the Himalayas and Hengduan Mountains have a relatively small 
pollen size (Table 1). Plants growing at high elevations with 
extreme conditions often have reduced morphological features 
[61,62], and pollen size also appears to show this tendency 
[33,63]. Furthermore, species native to Northeast China, such as 
L. leichtlinii var. maximowiczii, L. concolor var. pulchellum and L. 
dauricum, also have a relatively small pollen size (Table 1). A long 
winter and short growing season forces rapid completion of plant 
growth cycles. Hence, small pollen size also appears to be more 
efficient in tolerating relatively extreme conditions. In contrast, L. 
brownii and L. formosanum have relatively bigger pollen grains, likely 
because they occur in an environment suitable for slow 
development that allow pollen grains to grow to a larger size, as 
discussed by Muller [32]. Therefore, pollen size could be related to 
the selective pressure to adapt to environmental conditions. In 
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Section/ 
Subsection 



58/50 



85/80, 




L.davidii 
92/9a — L.pumilum 

' L.davidii var. willmottiae 

— L. concolor var. pulchellum 
L.leichtlinii var. maximowiczii 
L. lancifolium 
.dauricum 
L.cernuum 
— L.amabileAF088l97 

L.formosanum 
L. brownii 
L.taliense 
L.wardii 

L.sempervivoideum HM045467 
L.amoenum 
L.bakerianum 

L.tsingtauense 
j — L.martagon AF088203 
' — L.distichum 

L.martagon var. pilosiusculum 
L.lophophorum 
L.nanum AF088210 
L.henricii AB035285 



L.canademe AB020457 



L.columbianum AF090963 
100/100 I L.carniolicum EF042782 
L.bosniacum EF042789 

L.kesselringianum AM292429 
L.duchartrei AB035283 

L.speciosum var. gloriosoides FFM045461 



L.leucanthum var. centifolium 
■ L.sargentiae AB020443 
■L.leucanthum 
L.sulphureum 
L.regale 

Northolirion bulbuliferum HM045476 
Cardiocrinum gigantum 



0.02 




Figure 8. Mapping pollen exine sculpture characters on the ML tree of Lilium based on nuclear internal transcribed spacer (ITS) 
sequence data (adapted from Du et al. unpublished data). Values along branches represent bootstrap (BS) of ML and MP, respectively. A, B, C, 
D, Representatives of subsection 5a&5b. (A) L.pumilum (Jia 09001 ). (B) L.davidii (Du 1 1 008). (C) L. leichtlinii var. maximowiczii (Du 11010). (D) L. concolor 
var. pulchellum (Du 11012). E, F, Representatives of subsection Leucolirion 6b. (E) L.brownii (Wang 11006). (F) LJormosanum (Du 12008). G, 
Representatives of subsection 5c (G) LAaliense (Jia 09002). H, I, Representatives of subsection Lophophorum III. (H) L.sempervivoideum (Liang and 
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Zhang 1985). (I) Lbakerianum (Du 12003). J, K, Section Martagon. (J) L.distichum (Jia 07001). (K) LAsingtauense (Du 11006). L, Representatives of 
subsection Lophophorum I. (L) L. Lophophorum (Liang and Zhang 1985). M, Representatives of subsection Lophophorum II. (M) L. henricii (Liang and 
Zhang 1985). N, Representatives of section Pseudolirium. (N) Lcanadense (Kosenko, 1999). O, Representatives of section Liriotypus. (0) L. bosniacum 
(Muratovic et al. 2010). P, Representatives of section Sinomartagon I. (P) Lduchartrei (Wang 1 1005). Q, Representatives of section Archelirion. (Q) L. 
speciosum var. gloriosoides (Wu et al. 2007). R, Representatives of subsection Leucolirion 6b. (R) L. regale (Du 11015). 
doi:1 0.1 371 /journal.pone.0087841 .g008 



addition, smaller pollen grains appear to survive better in habitats 
with extreme conditions. 

Within the same species, pollen parameters from different 
provenances showed some variation (Table 1). Based on one-way 
ANOVA and LSD tests in 1 1 species (Table 2), results showed 
significant or extremely significant differences in pollen parameters 
in the same species from different habitats within the provenances, 
except for several single parameters in some species (Table 2). This 
may be explained as a consequence of interactions between 
phylogenetic and environmental constraints, as emphasized by 
Kawano and Kato [31]. 

Due to the wide variety of habitats in the distribution of Lilium, 
correlation analysis was performed to test which environmental 
factor is critical in influencing pollen parameters. The correlation 
analysis showed a significantly positive correlation between the 
polar axis and annual precipitation (P<0.05) and an extremely 
significantly positive correlation between the P/E ratio and annual 
precipitation (P<0.01) (Table 4). Hence, pollen size and shape 
showed a significantly positive correlation with annual precipita- 
tion. The absolute size of pollen grains influences harmomegathic 
function [64], which can be reduced to such a degree that exine 
elasticity alone can accommodate the necessary changes in volume 
to adapt to environmental conditions [30]. 

Pollen exine sculptures of different species in similar environ- 
mental conditions may be different types. For instance, L. concolor 
var. pulchellum possesses Concolor type pollen, whereas L. dauricum 
shows Martagon type pollen, and both occur in Northeast China 
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under similar environmental conditions. Species native to the 
Himalayas and Hengduan Mountains show all of the pollen types 
within Lilium. The evolutionary trend of exine sculptures is likely 
not definitively correlated with pollen size. 

Conclusion 

In conclusion, we suggest recognising a new pollen type, 
Formosanum, to accommodate pollen from L. formosanum. Pollen 
sculpture patterns appear to reflect phylogenetic relationships and 
are useful for species or subsection delimitation within sections. 
The reduction and simplification evolutionary trend of pollen 
sculpture and size could be related to the selective pressure to 
better adapt to environment conditions, especially extreme 
environmental conditions. In addition, the evolutionary trends of 
exine sculpture and pollen size and shape are not definitively 
correlated. Pollen size and shape show a significantly positive 
correlation with annual precipitation. However, additional studies 
are required to confirm these evolution hypotheses in a broad 
sample of the genus using techniques such as transmission electron 
microscopy (TEM) to observe strata of the exine. 
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